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Abstract 

Interest in oblique-wing aircraft has surfaced 
periodically since the 1940's. This concept 
offers some substantial aerodynamic performance 
advantages but also has significant aerodynamic 
and inertial cross-coupling between the aircraft 
longitudinal and lateral-directional axes. This 
paper presents a technique for synthesizing a 
decoupling controller while providing the desired 
stability augmentation. 

The proposed synthesis procedure uses the 
concept of a real model-following control sys- 
tem. Feedforward gains are selected on the 
assumption that perfect model-following condi- 
tions are satisfied. The feedback gains are 
obtained by using eigensystem assignment, and 
the aircraft is stabilized by using partial state 
feedback. The effectiveness of the control laws 
developed In achieving the desired decoupling is 
illustrated by application to linearized equations 
of motion of an oblique-wing aircraft for a given 
flight condition. 


Nomenclature 

A,B,C system matrices 

d eigenvector of components, unspecified 

e error 

I identity matrix 

j vn 

K feedback gain 

L,M matrices of appropriate dimension 

OWRA oblique-wing research aircraft 

p roll rate, deg/sec 

q pitch rate, deg/sec 

RMF real model -fol lowing 

r yaw rate, deg/sec 

♦Associate' Professor, Electrical Engineering 
Department. 

♦♦Aerospace Engineer. 


s 

complex frequency 

u ’ 

input vector 

V 

specified components of eigenvector 

w 

vector m-dimenslonal 

X 

state vector 

y 

output vector 

1 

eigenvector achievable for specified 
components 

a 

angle of attack, deg 

6 

sideslip angle, deg 

5 

control surface deflection 

e 

pitch angle, deg 


ei genvalue 

* 

bank angle, deg 

(...) R 1 

reordering operation 

Subscripts 


aL 

left aileron 

aR 

right aileron 

hL 

left horizontal 

hR 

right horizontal 

i 

i^h value 

m 

model 

P 

plant (aircraft) 

R 

set of real numbers 

u 

input vector 

X 

state vector 

Superscripts 

d 

desired value 

i 

number of outputs 



m number of inputs 

n number of states 

t pseudo-inverse 

Introduction 

The advantages of an oblique wing were first 
noted in the 1940's. However, only recently has 
the interest, technology, and mission of an 
oblique-wing design evolved into a full-scale 
flight research program. The NASA Ames Research 
Center and the U.S. Navy are developing an 
oblique-wing research aircraft (OWRA). Gregory* 
has outlined a number of potential advantages and 
disadvantages of this type of airplane. Theo- 
retical and wind-tunnel studies have shown that a 
variable skew oblique wing offers a substantial 
aerodynamic performance advantage for aircraft 
missions that require both high efficiency in 
subsonic flight and supersonic dash or cruise. 

The most obvious disadvantage of the oblique-wing 
concept is the asymmetry associated with wing-skew 
angle. This asymmetry results in significant 
aerodynamic and inertial cross-coupling between 
the aircraft longitudinal and lateral -di rectional 
axes. Current typical design procedures synthe- 
size aircraft controllers based on 2- or at most 
3-degree-of-freedom solutions. However, the OWRA 
stabilization and decoupling must consider at 
least 5 degrees of freedom simultaneously. 

The basis for OWRA will be NASA's F-8 digital 
fly-by-wire aircraft. This aircraft will be 
modified by the removal of the current high wing 
and installation of a wing pivot and a composite 
wing. A major part of the OWRA program will be 
the synthesis of a control system that will pro- 
vide acceptable stabilization and decoupling 
across the Mach, angle-of-attack , and wing-skew 
envelope. The aircraft thus offers an opportunity 
to apply modern control theory techniques to the 
solution of problems associated with OWRA. 


In this paper, the development of control laws 
for OWRA by Integration of the two above-mentioned 
techniques is demonstrated. The results show the 
effectiveness of the controller in obtaining the 
decoupled response for a given flight condition 
and wing skew. 

Problem Definition 

The concept of model -fol lowi ng is useful when 
an ideal set of plant equations of motion can be 
specified. The ideal objective of model -fol lowi ng 
flight control is to force the aircraft to respond 
as the model would to a given pilot command. It 
is often desirable to simulate the model dynamics 
in the flight computer and to generate the air- 
craft control signal using the aircraft outputs, 
the pilot input commands, and the model states. 
This situation is sometimes referred to as the 
pilot flying the computer, while the computer is 
flying the aircraft. 

More precisely, the model -fol lowing problem 
can be stated as follows. The linearized dynamics 
are given as 

x p = Ap x p + BpUp (1) 

Yp = CXp (2) 

where Xp e R", e R m , and y p e R A , and A p , Bp, 
and C are matrices of appropriate dimensions. The 
control Up must be determined such that the plant 
output yp approximates, reasonably well, some 
model output vector y m defined by the equations: 


x m " Am x m + B m u m 

(3) 

Ym = C x m 

(4) 


where x m e R n , u m e R m , and y m e R*, and A m , B m , 
and C are matrices of appropriate dimensions. 


Model -fol lowing has been a popular method for 
the design of multivariable control systems over 
the last two decades. In this method, the desired 
behavior of the plant is provided by an ideal 
model, and the problem is one of designing a 
suitable controller for the plant so that its 
response follows that of the model. 

Yore2 indicated the use of this method for 
simultaneous stability augmentation and mode 
decoupling. His synthesis procedure consisted of 
constructing an ideal model, designing feedback 
gains by quadratic optimization, and designing 
feedforward gains. A disadvantage of this method 
is that selection of feedback gain is an iterative 
and time-consuming process. The determination of 
this gain becomes a more complex problem when all 
states are not available and therefore output 
feedback is used. 

Another technique for decoupled flight control 
design is the eigenstructure assignment. 3 In 
this technique, the performance specifications 
can be interpreted in terms of the eigenvalues 
and eigenvectors of the closed-loop system. 
Broussard and Berry^ have established the 
equivalence of this technique to the design 
using model-following systems. 


For OWRA, the state, input, and output vectors 
are given by 



velocity, m/sec 
angle of attack, deg 
sideslip angle, deg 
bank angle, deg 
pitch angle, deg 
roll rate, deg/sec 
pitch rate, deg/sec 
yaw rate, deg/sec 


u 


{ hL 
fi hR 
«aL 
5 aR 
_«R . 


left horizontal tail deflection, deg 
right horizontal tail deflection, deg 
left aileron deflection, dpg 
right aileron deflection, deg 
rudder deflection, deg 


y 



0 


roll rate, deg/sec 
pitch rate, deg/sec 
yaw rate, deg/sec 
bank angle, deg 
angle of attack, deg 
sideslip angle, deg 
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The desired model of the aircraft, defined by 
matrices A m and B m as well as the aircraft matri- 
ces Ap and Bp are given in Table 1. The aircraft 
matrices correspond to a flight condition of 0.8 
Mach number and an altitude of 6096 m at 45° wing 
skew. The model used in this study is a modifica- 
tion of the zero-wi ng-skew configuration at the 
same flight condition. A m and B m elements are 

modified to increase damping and to eliminate 
zero-wi ng-skew coupling terms. This model is pre- 
liminary and may not represent ideal dynamics but 
does incorporate the desired aircraft decoupling. 

Mode 1 - F o 1 1 o w i n g Con t ro 1 System 

There are two configurations of model - 
following, one is implicit model -fol 1 owi ng , and 
the other is real model -followi ng (RMF). In 
implicit model-following, the model is not part of 
the system. In RMF, however, the model is part of 
the system as control law requires the states of 
the model. The technique of RMF has been shown to 
be amenable to the solution of many aircraft 
control problems. 

Frzberger 6 established conditions for perfect 
model -fol 1 owi ng that enable an ideal match of the 
dynamics of the compensated plant with those of 
the model. However, the conditions for perfect 
model -fol 1 owi ng are never attainable in the real 
world. An asymptotic RMF control law was derived 

by Chan 6 for the class of plants and models whose 
output vectors are identically their state vec- 
tors. Chan showed that, even if the conditions 
for perfect model -fol 1 owi ng are not satisfied, use 
of perfect model -fol lowi ng gains can yield a 
control capable of keeping error between the model 
and plant to a "small" region of statp space. 

Chan chose u p as 

Up = ui + U2 (5) 

where 

ui = K e (6) 

u 2 = Bp (A m - Ap)X|r| + BpB m u n 

= ^xm x m + Kum u m (7) 

and bJ is the pseudo-inverse of B p , and K xm and 
K um are the feedforward gains using model states 
and command input. Also, 

e = x m - x p (8) 

The control u p will ensure perfect model - 
following, if it is possible. If perfect model- 
following is not possible, the error setting 
rates would depend on eigenvalues of the closed- 
loop system and can thus be controlled in RMF. 
Also, if only partial state feedback is possible 
in the plant, perfect real model -fol 1 owi ng is 
still possible. 

For OWRA, because partial state feedback is 
to be used, the feedback gain K must be selected 
to ensure stability of the closed-loop aircraft 
and placement of its closed-loop eigenvalues 
at the desired location in the s-plane. The 


method of ei genstructure assignment will be 
used to select the gain K. This enables the 
desired eigenvectors and eigenvalues to be 
selected to ensure satisfactory tramsient 
response of the aircraft. 

Ejgensystem Synthesis 

Two widely used synthesis techniques of modern 
control theory are the linear quadratic regulator 
design and the modal control theory involving pole 
placement or eigenvalue-eigenvector assignment. ? 
One of the purposes of feedback control of air- 
craft is to improve or enhance the flying quali- 
ties of an aircraft. The difficulty in incor- 
porating specifications such as damping, natural 
frequency, and decoupling within a quadratic per- 
formance index makes the eigensystem synthesis 
procedure a promising design alternative. The 
performance specifications can be interpreted 
in terms of desired closed-loop eigenvalues and 
eigenvectors. Moore fi and others have shown how 
feedback can be used to place closed-loop eigen- 
value and shape closed-loop eigenvectors. Cun- 
ningham^, and Andry, Shapiro, and Chung, 10 and 

Sobel and Shapiro 11 have successfully demon- 
strated the use of ei genstructure assignment 
procedure for aircraft control system design. 

The handling qualities data base may be used 
to obtain desired pole locations directly. The 
additional design objective of obtaining augmented 
dynamics similar to those obtained in flight leads 
to specifications on eigenvectors or desired mode 
shapps. For example, pitch attitude must be domi- 
nant for the short-period mode, and speed must be 
dominant for the phugoid mode. 

Detailed discussions on eigenspace may be 
found in Ref. 11. However, some basic results for 
controllable and observable systems are summarized 
in the following discussion. 12 

Consider the system 

x = Ax + Bu 
y = Cx 

where x e R n , u e R m , and y e R*, and A, B, and C 
are matrices of appropriate dimensions. If the 
system is controllable and observable, and the 
matrices B and C are full rank, the following 
results hold 

1. The positions of maximum (m,£) closed-loop 

eigenvalues can be assigned arbitrarily with the 
stipulation that if X is a complex closed-loop 
eigenvalue, its complex conjugate X* must also be 
a closed-loop eigenvalue. 1 

2. The shape of maximum (m[i£)) eigenvectors 
can be altered. If the shape of a complex eigen- 
vector v-, is altered, its complex conjugate v* 
must be altered in the same way. 

3. For each eigenvector whose shape is 
altered, minimum (m,r) eigenvector elements 
can be chosen arbitrarily. 
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4. Attainable eigenvectors must lie in the 
subspace spanned by the columns of ( X-, I - A)"*B 
of dimension m that is the number of independent 
control variables. A desired eigenvector v^ will, 

in general, not reside in the prescribed subspace 
and cannot be achieved. The achievable eigen- 
vector v* is obtained by orthogonal projection 

of vj onto the subspace spanned by ( X-, I - A)’^B. 

It will generally be true that only a few of the 
components in v 1 are actually specified. The 
remainder can be arbitrary. To account for this, 
v ] is reordered and partitioned as follows 

. (£-) (9) 

where 


v* = the specified subvector 

d-, = the vector of unspecified components 
(. . .) R i * the reordering operation 

If we let 

(£-)= [<m - ArVK do 

then, as shown in Ref. 13, Z-, may be selected to 
best approximate v* with v°. By the method of 
orthogonal projections, Z-, is obtained 

z, = (L'L)" 1 L'v* (11) 

As shown by Moored, the feedback gain K is 
given by 

K = (w 1 w 2 ...w n )(viV2...v n )' 1 (12) 

where w-| is obtained from the relation 

(At I - A) v-| = Bw n (13) 

Results 

To illustrate the degree of coupling in the 
open-loop system and decoupling in the closed- 
loop system, a one-degree control command was 
input for 2 sec as shown in Fig. 1. This command 
input was either elevator or aileron and was 
reduced to zero after 2 sec. Figures 2(a) to 2(c) 
illustrate the open-loop system response to an 
elevator command input for pitch rate, yaw rate, 
and bank angle, respectively. Significant yaw 
rate and bank angle are generated as a result of 
the pitch command, and of particular interest is 
the very large change in bank angle illustrating 
the significant cross-coupling. 

Table 2 shows the desired eigenvector assign- 
ment specification, open- and closed-loop eigen- 


values and desired eigenvalues, the feedback 
gain matrix K, and the feedforward gain matrices 
K xm and K um . 

Figures 3(a) to 3(c) and 4 illustrate the 
closed-loop system response to the same elevator 
command input. The pitch rate in Fig. 3(a) is 
attenuated as compared with the open-loop 
response. However, the system is very closely 
following the desired model response as illus- 
trated in Fig. 4. The yaw rate and bank angle are 
virtually nonexistent as illustrated in Figs. 3(b) 
and 3(c), thus achieving the desired decoupling. 

Figures 5(a) to 5(c) illustrate the open-loop 
system response to a one-degree aileron command 
input. Pitch and yaw angular rate and bank angle 
are again shown. The relative pitch coupling is 
not as severe for the aileron command case as the 
roll coupling is for the elevator command, how- 
ever, coupling is still present. 

Figures 6(a) to 6(c), 7(a), and 7(b) illus- 
trate the closed-loop system response to the 
same aileron command. Pitch rate is virtually 
nonexistent, and the desired yaw rate and bank 
angle are achieved, giving the desired decoupling. 

Conclusions 

A method is presented to obtain a decoupled 
control for a highly coupled asymmetric aircraft. 
The method utilizes a real model-following control 
law in which gams for perfect model-following are 
used even when the conditions for perfect model - 
following are not satisfied. The feedback gain, 
using output feedback, is computed by using 
eigenstructure assignment. The results indicate 
that the method does obtain the decoupling incor- 
porated in the ideal model for the flight con- 
dition considered. 

Future investigations will be conducted to 
evaluate the control algorithm under nonlinear 
6-degree-of-f reedom flight conditions. These 
investigations will consider such factors as 
nonlinear aerodynamic data, control system 
surface rate and position constraints, and 
system hysteresis. 
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TABLE 1. - AIRCRAFT AND MODEL MATRICES 





Plant 

(ai rcraft) 

matrices 





r 0.0094 

22.0707 

10.5479 

-0.1341 

-32.1127 

0.0057 

-0.0005 

-0.0265" 


-0.0001 

-0.7826 

0.0958 

0.0000 

0.0000 

0.0030 

0.9926 

-0.0003 


0.0000 

-0.0592 

-0.2908 

0.0387 

-0.0002 

0.0259 

0.0001 

-0.9920 

A p = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0247 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

33.1432 

-53.6933 

0.0000 

0.0000 

-3.1250 

2.0552 

1.7210 


-0.0002 

-8.6816 

0.7975 

0.0000 

0.0000 

0.1679 

-1.0352 

0.1810 


L o.oooo 

-1.0092 

10.7521 

0.0000 

0.0000 

-0.0213 

0.0080 

-0.7129 J 


r 0.0844 

-0.0309 

-0.2210 

-1.2572 

3.65981 





-0.0974 

-0.0974 

-0.0198 

-0.0302 

0.0000 





-0.0166 

0.0166 

0.0008 

-0.0005 

0.0647 




Bp = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 




0.0000 

0.0000 

0.0000 

0.0000 

0.0000 





12.9804 

-22.2654 

15.8467 

-11.8422 

13.2774 





-9.4073 

-10.8655 

-1.2311 

0.8797 

0.5694 





L 1.9854 

-2.2579 

0.5262 

-0.3276 

-6.2499^ 




Model matrices 


r-0.0077 

23.5966 

0.0000 

0.0000 

-32.1129 

0.0000 

0.0000 

o.ooool 


-0.0001 

-1.1062 

0.0000 

0.0000 

0.0000 

0.0000 

0.99U9 

0.0000 


0.0000 

0.0000 

-0.6000 

0.0387 

0.0000 

-0.0148 

0.0000 

-0.9919 

A 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

-0.0133 

A m - 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

0.0000 

-44.3777 

0.0000 

0.0000 

-10.0000 

0.0000 

0.0000 


0.0001 

-12.1514 

0.0000 

0.0000 

0.0000 

0.0000 

- .0000 

0.0000 


L o.oooo 

0.0000 

12.1943 

0.0000 

0.0000 

0.0000 

0.0000 

-2.0000 ij 


[“ -2.2032 

-2.2032 

-0.8354 

-0.8354 

0.0000 1 





-0.0848 

-0.0848 

-0.0494 

-0.0494 

0.0000 





-0.0166 

0.0166 

0.0000 

0.0000 

0.0647 




Bm = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 




u m 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 





11.2379 

-11.2379 

29.0513 

-29.0513 

9.6847 





-7.8229 

-7.8229 

0.0000 

0.0000 

0.0000 





_ 0.0000 

0.0000 

0.0000 

0.0000 

-6.6502J 
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TABLE 2. - EIGENSYSTEM ASSIGNMENT AND GAINS 


Desired eigenvectors 3 


Short period Dutch roll 



Spiral Roll subsidence 


"0“ 


"0” 

0 


0 

0 


0 

1 


X 

0 


0 

X 


1 

0 


X 

-X. 


_0_ 


Eigenvalues 


Condition 


Open loop 


Desired Achieved 

closed loop closed loop 


Short period 
Dutch roll 
Spi ral 

Rol 1 subsidences 
Phugoid 


-1.0433 ± j 2 • 8269 
-0.5463 ± j 3 . 381 6 
-0.0118 
-2.7544 

-0.0053 ± jO. 0455 


-2 ± j 3. 5 -2 ± j 3.5 

-3 ± j4.0 -3 ± j 4.0 

- 0.1 - 0.1 

-7.0 -7.0 

-0.0047 t jO.0455 


Feedback gain 



p 

q 

r 

* 

a 

0 



” 0.1454 

-0.3471 

-1.4507 

0.1476 

-0.5902 

-0.95871 



-0.0580 

0.4180 

1.0127 

-0.0915 

1.0377 

0.1151 


K = 

-0.3710 

0.9639 

2.2955 

-0.2967 

-5.0019 

6.9733 



0.1199 

0.1133 

-0.0982 

0.0194 

-7.2635 

1.3083 



L 0.0114 

-0.1653 

0.1540 

0.0210 

0.6898 

-1.9884J 




Feed forward 

gams 






K 

xm “ Bp(Am 

■ Ap ) 




V 

a 

0 

4> 

9 

P 

q 

r 


r o.ooio 

6.5743 

-3.6200 

0.0729 

-0.0003 

0.1698 

-0.2893 

-0.359f 


-0.0006 

-3.6600 

2.2241 

-0.0452 

0.0002 

-0.0560 

0.4205 

0.2696 

K xm " 

-0.0019 

-9.0867 

9.4682 

-0.1429 

0.0004 

-0.5446 

0.3861 

0.5095 

0.0002 

7.0692 

2.0716 

0.0037 

-0.0002 

0.1608 

-0.6155 

-0.0284 


L 0.0004 

2.1139 

-1.4962 

0.0272 

-0.0001 

0.0164 

-0.1778 

0. 0388_ 



^urn 

= BpB m 



”-0.8204 

-0.5391 

-0.0271 

0.7842 

-2.1203“ 

1.0516 

0.9978 

-0.1411 

-0.3323 

1.3850 

3.1869 

1.6393 

2.1149 

-1.9684 

3.9478 

0.0239 

0.2235 

0.8121 

1.4435 

-0.1933 

.-0.3735 

-0.4054 

0.1779 

0.1278 

0.2326 


a x is arbitrary. 


Elevator 

deflection, 

deg 



Fig . 1 Command input to system . 
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(a) Pitoh rate. 


(b) Yaw rate , 



(o) Bank angle , 

Fig , 2 Open-loop system response to elevator oormand input , 



(aJ Pitoh rate , 


Yaw 

rate, 

deg/sec 




fW yaw rate. Bank angle. 

Fig, 3 Closed-loop system response to elevator command input. 
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Time, sec 
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Fig . 4 Model-following response to elevator 
command input • 




4 5 l 

Time, sec 

Yaw rate . 


(c) Bank angle. 

Fig * 5 Open-loop system response to aileron command input 



Pitch 

rate, 

deg/sec 


0 


01 23456789 10 

Time, sec 

(a) Pitch rate . 



(b) Yaw rate . (c) Bank angle . 


Fig • Closed-loop system response to aileron command input . 


6 
5 
4 

Yaw 3 
rate, 

deg/sec 2 
1 
0 

- 1 

0123456789 10 

Time, sec 

(a) Yaw rate . 


Plant (aircraft) 

Model 



Bank angle . 


Plant (aircraft) 

Model 



Fig. 7 Model-following response to aileron command input . 
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